function prob=p_bounded(x,beta,c,f)
	prob = sigma(f)+(sigma(c)-sigma(f))*(1/(1+exp(-dot(x,beta))));   
end
